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Principal  Nomenclature 

A1 

a2 

a  -  (/$“  -  J §y)l2 
b  = 

C 


Cp(0) 


ds 


E 

F  **  <J>  +  iijj 


’I' 


C 

O' 


U 


w 

,  dF 
w  =  u  -  iv  =  t— 
dz 


Z 


cavity  detachment  point  near  or  at  profile  nose 

cavity  detachment  point  at  profile  trailing  edge 

see  related  nomenclature  below 

see  related  nomenclature  below 

location  of  eigensolution  singularity  on  unit 

circle 

drag  coefficient 
lift  coefficient 

pressure  coefficient  on  wetted  surface 

profile  chord  length,  c  =  1 

element  of  arc  length  in  z-plane 

strength  of  eigensolution 

complex  potential  and  complex  F-plane 

velocity  potential 

stream  function 

stagnation  point  location 

point  at  infinity  in  z-plane 

free-stream  velocity 

intermediate  mapping  complex  plane 

complex  velocity 


z 


X  +  iy 


intermediate  mapping  complex  plane 
complex  variable  in  the  physical  (x,y)  plane 
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C 

n 

£ 

o 

T 
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Yc 


*1 

<f>2 

oi(c) 


angle  of  attack 

angular  location  of  stagnation  point  measured 
from  negative  real  axis 
dummy  variable 

normal  effect  from  profile  chord 

complex  semi-circle  plane,  £  =  e*-® 

distance  along  profile  chord 

cavity  thickness  at  trailing  edge 

stagnation  point  angular  location  on  unit 

circle  in  C  plane  (cos  y  =  a/b) 

angular  location  of  eigensolution  singularity 

on  unit  circle 

value  of  at  Aj 

value  of  $  at  A2 

9  +  i£n  q/U  =  9  +  it  complex  logarithmic  hodograph 

9  flow  inclination 

q  flow  speed 


Subscripts 

0 

1 


c 


pertains  to  flat  plate  solution  ujq 
pertains  to  regular  part  of  solution  wj 
pertains  to  eigensolution  toc 
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Introduction 

The  present  paper  bears  upon  the  two-dimensional  inverse  or  design 
problem  for  fully-cavitating  hydrofoils  in  which  one  specifies  the  pressure 
distribution  on  the  profile  wetted  surface  and  then  calculates  that  wetted 
surface  shape  which  will  satisfy  this  prescription.  This  design  problem 
is  certainly  not  new  to  airfoil  designers  and  as  far  as  cavity  flows  are 
concerned,  both  linear  and  nonlinear  design  methods  have  been  worked  out. 

In  the  realm  of  nonlinear  approaches  to  the  problem,  the  very  general 
method  of  Yim  and  Higgins  [1]*  is  worthy  of  note  because  it  applies  to 
single  foils  as  well  as  to  cascades  of  profiles  for  all  cavitation  numbers 
in  the  cavity-flow  regime.  Another  approach  has  been  discussed 
superfically  by  Khrabov  [2].  Both  of  these  contain  far  more  generality 
than  is  required  for  this  study  at  zero  cavitation  number.  For  the 
direct  or  off-design  problem  of  exact  cavity-flow  theory,  a  good  example 
of  the  present  level  of  development  is  represented  by  the  work  of  Furuya 
[3]  and  it  is  clear  that  now  one  can  do  both  the  design  and  off  design 
problems  for  fully-cavitating  hydrofoils.  Thus,  one  can  attempt  to 
tailor  the  profile  to  an  entire  set  of  performance  goals  and  failing 
that  he  can  at  least  design  for  the  best  compromise  among  a  set  of 
conflicting  requirements. 

According  to  many  authors  {4-7] ,  the  inverse  problem  is  not  thought 
to  present  much  of  a  challenge  at  zero  cavitation  number.  In  this  case, 
the  classical  method  of  Levi  Civita  [7]  can  be  applied  to  an  isolated 


•Jf 

Numbers  in  square  brackets  indicate  citations  in  the  references  listed 
below. 
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body.  This  view  is  certainly  proper  as  long  as  one  is  content,  after 
prescribing  the  pressure  in  the  circle  plane,  to  accept  whatever 
correlation  between  points  in  the  circle  and  physical  planes  may  result. 

Of  course,  such  a  rudimentary  approach  does  not  lead  to  a  useful  design 
procedure . 

The  motivation  for  the  present  investigation  is  that  none  of  the 
literature  on  the  nonlinearized  inverse  problem  we  have  surveyed  so  far 
[1-8]  has  made  use  of  the  fact  the  exact  theory  admits  the  existence  of 
a  "point  drag"  or  complementary  function.  While  this  fact  is  of  no 
particular  importance  in  the  direct  problem,  we  have  already  seen  in 
the  case  of  the  linearized  inverse  problem  [9-11]  that  the  complementary 
function  can  play  an  important  role.  For  the  exact  theory  there  has  been 
a  question  if  a  nonlinear  eigensolution  exists  [1].  Therefore,  in  this 
study,  we  explore  the  question  regarding  the  existence  of  a  "point  drag" 
or  eigensolution  in  the  nonlinear  theory  under  the  simplest  set  of 
circumstances  and  this  leads  us  naturally  to  the  restrictions  that  the 
free  streamline  flow  pertains  to  an  isolated  profile  and  that  the  flow 
be  at  zero  cavitation  number.  These  simplifications  will  free  us  from 
the  complications  arising  from  non-zero  cavitation  numbers  and  other 
boundaries  in  the  flow  domain  such  as  a  free  surface  or  neighboring 
cascade  blades. 

In  this  paper  we  use  the  term  eigensolution  in  the  sense  of  thin 
airfoil  theory  as  suggested  by  the  work  of  Van  Dyke*  because  we  already 
know  that  the  inverse  problem  in  the  theory  of  fully  cavitating  hydrofoils 

*Perturbatlon  Methods  In  Fluid  Mechanics,  The  Parabolic  Press,  Stanford, 
CA,  pp.  48-54  (1975). 
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is  not  necessarily  unique.  Our  aim  is  to  find  a  sufficiently  weak 
singularity  which  can  be  added  to  the  classical  Levi  Civita  solution  and 
can  then  be  used  to  satisfy  certain  additional  physical  conditions  relating 
to  the  location  the  free  streamline  springing  from  the  hydrofoil  nose 
and  thereby  provide  a  unique  inverse  cavity-flow  solution.  After  we  have 
constructed  the  simple  eigensolution,  we  will  examine  some  of  its 
properties.  The  use  of  this  solution  in  the  design  process  will  be 
presented  elsewhere. 

Flow  Geometry  and  Conformal  Mappings 

As  noted  above,  this  study  uses  Levi  Civita 's  method  [7]  for  the 
analysis  of  the  exact  inverse  problem  for  a  fully  cavitating  hydrofoil 
section.  The  flow  geometry  and  the  principal  quantities  associated  with 
the  flow  are  illustrated  in  Fig.  1.  The  origin  of  coordinates  in  the 
z  =  x  +  iy  plane  is  taken  at  the  stagnation  point  on  the  wetted  surface 
of  the  hydrofoil.  This  point  is  denoted  by  0  in  Fig.  1.  The  chordline 
of  the  profile  is  inclined  at  the  angle  a  with  respect  to  the  x  axis 
and  the  free-stream  velocity  U  is  taken  as  being  parallel  to  this  axis 
as  illustrated.  The  flow  separation  point  at  or  near  the  profile  nose 
is  the  point  Aj  as  illustrated  for  a  sharp-nosed  foil.  The  "upper’' 
cavity  surface  is  shown  as  the  dashed  curve  extending  from  A^  to  the 
point  0  at  downstream  infinity.  In  the  case  of  a  round-nose  profile 
Aj  can  lie  on  the  upper  wetted  surface  behind  the  leading  edge.  This 
case  is  not  illustrated  in  Fig.  1.  The  point  A2  denotes  the  location 
of  the  trailing  edge  of  the  wetted  surface.  The  lower  surface  of  the 
cavity  leaves  the  wetted  surface  at  A2  and  extends  as  shown  by  the 
dashed  line  to  the  point  O'  at  downstream  infinity. 
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Let  the  coordinates  of  a  typical  point  on  the  wetted  surface  be  denoted 
by  x  and  y  and  those  on  the  upper  cavity  by  xc  and  y^  as  shown  in  Fig.  1. 
While  the  orientation  of  the  profile  in  the  z  plane  is  convenient  for  the 
purposes  of  analysis,  the  x-y  system  is  not  always  a  convenient  reference 
frame  for  foil  and  cavity  contours.  For  this  purpose  we  use  a  coordinate 
system  with  the  abscissa  along  the  chordline  as  shown  by  the  distance  a 
measured  from  the  profile  nose.  The  ordinates  of  the  wetted  surface  are 
then  given  in  terms  of  a  as  n(a)  and  the  upper  cavity  ordinates  are  given 
by  nc(o).  At  the  trailing  edge  of  the  profile,  the  cavity  thickness  is 
nc  =  T.  These  quantities  are  also  shown  in  Fig.  1.  In  the  o,n  system 
the  stagnation  point  0  is  located  at  (ag»hg)  as  illustrated.  The  trans¬ 
formation  between  the  (x,y)  and  (a,n)  systems  is 

a  +  in  -  aQ  +  ihg  +  zeia  ,  (1) 

where  z  is  the  complex  variable,  z  «  x  +  iy,  and  a  is  the  angle-of-attack 
as  measured  by  the  inclination  of  the  chordline  with  respect  to  the  x  axis 
and  free-stream  velocity  U. 

The  conformal  mappings  start  with  the  complex  potential  in  the 
z-plane , 

F  -  $  +  U>  ,  (2) 

where  <j>  is  the  velocity  potential  and  ^  is  the  stream  function.  As  is 
customary,  we  adjust  these  quantities  to  make  $  =  0  at  the  stagnation 
point,  0.  The  stream  function  is  taken  to  be  zero  all  along  the 
stagnation  streamline.  Therefore,  the  boundaries  of  the  flow  can  be 
represented  by  a  cut  all  along  the  real  axis  in  the  F-plane  as  shown  in 
Fig.  2.  Note  that  the  wetted  surface  extends  from  the  stagnation  point 
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at  0  to  the  trailing  edge  at  A2  and  the  lower  cavity  surface  from  A2  to  o' 
must  lie  along  the  lower  surface  of  the  cut.  This  is  so  because  downstream 
from  the  stagnation  point,  the  velocity  potential  increases  in  the  flow 
direction  and  the  stream  function  decreases  outward  from  the  foil  or 
cavity  surface.  On  the  arc  OAj ,  the  flow  direction  is  reversed  with  a 
consequent  reversal  in  the  gradients  of  <J>  and  so  that  the  point  Aj  is 
on  the  upper  edge  of  the  cut. 

One  can  use  the  mapping, 

W  =  ✓?  (3) 

in  order  to  map  the  F-plane  outside  the  cut  into  the  upper  half  of  the 
W-plane.  Corresponding  points  are  shown  by  the  locations  of  o',  Aj  and 
A2  and  0  in  Fig.  3.  As  before,  the  cavity  surfaces  are  shown  as  dashed 
lines  and  the  wetted  surface  is  shown  by  the  solid  line.  Let  the  values 
of  4>  at  Aj  and  A2  in  the  F-plane  be  4>  1  and  $2  respectively.  Then  these 
points  are  at  W  =  /<jT|~  and  W  =  -  and  the  midpoint  of  the  distance 
between  Aj  and  A2  is  located  at 

W  =  a  =  — -i-j - -  .  (4) 

The  distance  between  this  midpoint  and  Aj  is 

/<t>,  + 

b  =  — - - -  .  (5) 

These  distances  are  also  shown  in  Fig.  3. 
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dF  iU)c  d(/F2)  1E//F 

az  ~  U  6  U  6 


If  we  put  t  =  E//F,  we  find  for  a  profile  of  unit  chord  that 


2  D 

,  2E  1  it.  c  1  it , 

dz  - - - - -  e  dt  - - re  dt 

"  t3  '  t3 


which  has  the  indefinite  integral, 


z 


(36) 


Completion  of  this  integration  can  be  carried  out  in  four  parts  starting  from 
either  side  of  the  stagnation  point  where  z  =  0  and  t  -*• 

For  example,  on  the  arc  OA^  we  recall  that  6  >  y.  Put  t“ti,Y*T-6 
and  0  j  =  tt  —  5 1  then 


t 


1 


_ E _ 

b(cos  -  cos  6) 


Because  0  <  C]  <  6,  we  have 

E 

- - -  <  t  <  “  . 

b(l  -  cos  o)  1 

After  separating  Eq .  (36)  into  real  and  imaginary  parts,  we  find  for  this 
part  of  the  wetted  surface,  the  coordinates 
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r  =  0 
c 


Since  tc  =  0  on  the  wetted  surface,  we  expect  ioc  to  make  no  contribution 
to  the  lift  although  the  singularity  at  the  stagnation  point  should  lead  to 
a  drag  force.  Making  use  of  Eqs.  (21),  (22)  and  (34a),  we  find  that 


Along  the  arc  OAj .  Therefore  9C  changes  sign  at  B  =  Y • 

The  shape  of  the  wetted  and  cavity  surfaces  follows  from  the 
relationships  of  Eq .  (14)  which  can  be  expressed  as 
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0) 

c 


0  +  It 

c  c 


E_ 

/F 


(33) 


where  E  is  a  real  constant,  we  have  a  function  which  satisfies  the  necessary 
requirements.  The  two  conditions  u)c(<*>)  =  0  and  qc  =  U  on  both  the  cavity  and 
the  profile  wetted  surface  can  be  satisfied  by  any  member  of  the  family  of 
functions  having  the  form  F_m,  0  <  m  <  1.  But  the  condition  0C  =  0  along  the 
entire  stagnation  streamline  can  be  satisfied  only  when  m  =  1/2.  This  choice 
for  the  complementary  function  seems  to  offer  the  advantage  that  it  will 
cause  less  alteration  of  the  upstream  flow  field  inclination  than  other 
possibilities.  Moreover,  it  appears  to  be  the  most  convenient  choice  for 
further  analysis.  Consequently,  we  shall  adopt  this  functional  form  for 
the  eigensolution  in  this  work. 

Equation  (10)  can  be  used  to  obtain  the  corresponding  representation 
for  <jc  in  the  5-plane  as 


»  =  _ l _  7  2EC 

b[cos  Y-y  (?+“■)]  b(c  "  e1Y)(c  -  e  1Y) 


(34a) 


In  the  5-plane  <uc(0)  =  0  and  when  5  is  real  mc  is  real.  Moreover,  u>c(5)  is 

an  analytic  function  which  is  regular  inside  the  unit  circle  and  which  has 

simple  poles  at  5  *  Note  from  Eq.  (34a)  that  on  the  "nose  cavity" 

0C  +  0+  as  5  +  -  0  and  on  the  "tail  cavity"  9C  ■+  0-  as  5  -*■  +  0,  as 

-iw 

illustrated  in  Fig.  1  also.  From  Eq.  (14),  w^  =  Ue  and  we  see  that 
the  structure  of  u>c  leads  to  an  essential  singularity  in  wc  at  the 
stagnation  point  0.  The  complex  velocity  wc  is  bounded  at  this  point 
however,  and  1  smooth  foil  contour  will  pass  through  z  =>  0.  At  points 
on  the  unit  circle,  5  =  e1-®,  we  have 
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These  results  will  be  used  in  order  to  start  the  functional  iteration  after 
the  formulation  of  the  theory  has  been  completed.  As  a  final  note,  we 
observe  that  the  cavity  thickness  at  the  trailing  edge  of  a  flat  plate 
according  to  linearized  theory  [13]  is 

nc ( 1 )  =  1.681a  , 

which  can  be  compared  to  the  corresponding  expression  from  Eq.  (31b).  Within 
the  range  of  attack  angles  considered  here  the  trailing  edge  cavity  thickness 
at  zero  cavitation  number  is  about  the  same  when  estimated  by  linearized 
or  nonlinear  theory.  For  larger  values  of  attack  angle  we  would  expect 
estimates  from  linearized  theory  to  exceed  nonlinear  cavity  thicknesses. 

The  flat-plate  function  wqCc)  is  traditionally  considered  to  possess 
all  of  the  singular  behavior  of  the  function  0>(O»  The  shape  of  the  smooth 
body  is  then  represented  by  an  analytic  function  mj(c)  which  is  regular 
inside  and  on  the  unit  circle.  It  must  also  satisfy  the  same  symmetry 
requirements  as  are  imposed  upon  wq  and  we  must  also  insist  that  u>j(0)  =*  0. 
Then  one  traditionally  puts  00(5)  =  ojq  +  a)j.  As  mentioned  previously,  we 
will  add  to  this  customary  sum  a  new  function,  cuc(c),  which  is  the  analog 
of  the  point-drag  function  of  linearized  theory.  We  will  now  explore  the 
properties  of  this  eigensolution. 

The  Eigensolutlon 

The  complementary  function  u>c  ■  9C  +  itc,  is  to  be  determined  from  the 
requirements  that  xc  *  0  on  the  cavity  and  the  foil,  9C  =  0  along  the 
stagnation  streamline  and  that  <uc  vanishes  at  infinity.  A  function  of  the 
form  which  satisfies  these  conditions  can  be  found  most  easily  by  considering 
the  flow  in  the  F-plane,  Fig.  2.  If  we  take 
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Equation  (28)  gives  the  value  of  to  be  used  in  Eqs.  (29)  and  (30). 

It  is  useful  to  find  the  ordinate  of  the  upper  surface  of  the  cavity 
above  the  trailing  edge  of  the  wetted  surface.  This  can  be  done  by 
calculating  the  quantity 


cos  a  - 


for  several  values  of  and  prescribed  values  of  a.  Then  when  ac  =  1 ,  one 
can  plot  contours  of  a  »  constant  in  an  f  -  £c  plane  and  note  for  each  value 
of  a  the  value  of  Cc  corresponding  to  f  =  2(4  +  n  sin  a).  These  values  of 
Cc  can  then  be  used  in  Eq.  (30)  in  order  to  compute  the  value  of  nc(l)  for 
each  value  of  a.  It  was  found  when  these  points  were  plotted  in  an 
nc(l)  -  a  plane  that  a  linear  relationship  fit?  the  data  for  0  <  a  <  10°. 
When  a  is  measured  in  degrees  this  line  has  the  equation 


nc(l)  =*  .2939a° 


(31a) 


The  corresponding  relationship  for  a  measured  in  radians  is 

nc(l)  =  1.684a  .  (31b) 


The  computed  results  are  tabulated  below. 


Table  1.  Cavity  Thickness  at  the  Trailing  Edge  of  a  Flat  Plate 


a° 

nc(l) 

0.0 

- . .  _  _ 

.0000 

0.5 

.1712 

.0147 

1 .0 

.1707 

.0297 

2.0 

.1698 

.0585 

3.0 

.1689 

.0879 

4.0 

.1680 

.1174 

5.0 

.1669 

.1463 

7.0 

.1650 

.2049 

10.1 

.1620 

.2939 
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Additional  properties  of  the  flat-plate  solution  which  are  important 
for  the  present  considerations  relate  to  the  shape  of  the  upper  surface  of 
the  cavity.  In  Fig.  (5)  5  will  be  on  the  negative  real  axis  for  points 
on  this  part  of  the  cavity.  Therefore,  let  us  put  £  *  -  where  ;c  is 
a  real  positive  number.  Then  we  can  use  d(o  +  in)  =  eiotdz,  and  write 
Eq.  (17)  in  the  form 

2 

d(o  +  in)  -  Jp  exp  i[o  +  fl0(5c)]d{[|  (?c  +  -  cos  a]2}  . 

But  now 

exp  i[a  +  eQ]  =  e"ia(eia  -  ?c)/(e_iCt  -  cj 

and  so  we  have 

2 

d<0  +  in)  -  y  [c  e"i0t  -  2  +  (21  sin  a)/?  +  Vl  -  eia/^3]dC  . 

z  u  c  c  c  c  c 

The  integration  in  this  case  starts  at  Aj  where  Cc  =  1  and  a  *  n  =  0  and 
proceeds  to  some  value  of  oc,nc  corresponding  to  1  <  J  <  0.  This  leads  to 
the  parametric  representation  of  nc(oc)  in  terms  of  ?c  which  is  given  by 


o 

c 


1 

2  U 


(29) 


and 


n 


c 


1  b 

2  U 


t  rr~  sin  a 


+  2  i  nc 

c 


(30) 


9 
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2 

,  b  (cos  £  -  cos  a)  i9rr  ,  \2  .  e  ,  ,  .  \2  ,  ,  ,, 

dz  -  -  - — [ - - - r  e  I  cos  £  ~  cos  a  +  I  sin  £  +  sin  a  sin  £  d£ 

U  cos  £  -  cos  a  v  ' 


i  0  b**  2 

*  -  e  sgn(cos  £  -  cos  a)  2  —  [sin  £  +  sin  a  sin'1  £  -  cos  a  cos  £  sin  £]d£ 

We  can  dispose  of  the  product  -  e^-®  sgn  (cos  £  -  cos  a)  by  observing  that  on 
the  wetted  surface  when  0  <  £  <■  a  the  flow  direction  is  0  =  ir  -  a  and  when 
a  <  £  <  ir,  9  =  -  a.  Therefore  in  the  first  of  these  cases,  e1-®  =  -  e“ia  and 
in  the  second  e^®  *  +  e“^a.  Hence  in  either  of  them  the  product 
-ei9  sgn(cos  £  -  cos  a)  =  e-^-01.  Next  we  can  introduce  d(o  +  in)  =  e^adz 
from  Eq.  (1)  with  the  result  that 

b2  2 

d(a  +  in)  *  2  —  [sin  £  +  sin  a  sin  £  -  cos  a  cos  £  sin  £]d£ 

This  last  result  implies  that  n  3  0  as  is  proper  for  a  flat  plate  and 
da  is  simply  the  arc  length  ds  along  the  wetted  surface  measured  from  the 
profile  nose  where  a  =  0.  We  can  integrate  this  last  equation  from  a  ■  0 
(£  =  0)  to  some  value,  0  <  o  <  1  (0  <  £  <  tt),  and  get 

s  -  2  ~  [(l  -  cos  £)(l  -  [i  +  cos  5])  +  (g  -  Sin2  g)  Si^_SL]  .  (27) 


However,  the  profile  is  to  have  unit  chord  so  that  when  £  =  n,  a=s=l  and 


U  4  +  it  sin  a 


(28) 


When  Eq.  (28)  is  used  in  Eqs.  (24)  and  (25),  one  obtains  the  well  known 
Rayleigh  formulae  [12]. 


-15- 


25  January  1985 
BRP : lhz 


2  cos  8  -  cos  y 

2  2 
(cos  y  -  cos  8)  +  (sin  y  +  sin  8) 


Next  we  can  put  8  =  it  -  £  and  y  =  n  -  a  in  this  expression  and  in  Eq.  (17) 
when  it  is  written  for  points  on  the  wetted  surface.  That  is, 
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the  real  part  of  the  analytic  function  i£n(c  -  e*Y).  However  this  function 
by  itself  does  not  have  the  proper  symmetry  needed  for  admissible  forms  of  to 
But  if  we  subtract  from  it  a  similar  function  which  has  a  like  jump  at  the 
image  of  0  with  respect  to  the  real  axis  in  the  ? -plane,  we  preserve  the 
necessary  behavior  at  0  and  also  satisfy  the  preceding  symmetry  requirements 
Aside  from  arbitrary  additive  constants,  this  function  has  the  form 


i£n 


Finally,  we  require  that  <o(0)  =  0.  Because  of  this  condition  the  resulting 
function  which  provides  the  flat  plate  solution  is  [4] 


WqU)  -  0Q  +  iTQ  =  -  a  +  Un 


C  -  e 


iY 


1  -  te 


iY 


(23) 


with  y  =*  it  -  a  in  the  case  of  an  isolated  flat-plate  profile.  For  this  case 
one  can  show  that  this  function  has  the  flow  direction,  9q  *  -  a  or  0q  *  ir  - 
on  C  *  e1-^.  Therefore,  the  wetted  surface  in  the  z-plane  is  a  straight  line 
through  0  with  its  trailing  edge  inclined  at  the  angle  -  a  with  respect  to 
the  positive  real  axis.  Moreover,  since  Tq  vanishes  on  the  real  axis  in  the 
C-plane,  the  free  streamlines  have  Cp  *  0  as  required. 

The  contribution  of  dig  to  Cp  and  Cp  follows  from  Eqs.  (21)  and  (22). 
From  Eqs.  (21),  (22)  and  (23)  we  get 

b2 

CT  =  2tt  — —  sin  a  cos  a  (24) 
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t(-8)  -  -  t(B) 

and 

0(-B)  -  0(8)  . 

Hence  t  is  an  odd  function  of  8,  or  Im  C,  and  0  is  an  even  function  of  8,  or 

Im  £. 

Using  the  result. 


eT  -  /I  -  C  (20) 

P 

as  noted  above,  and  the  fact  that  u)(0  is  now  defined  inside  the  unit  circle 
one  finds  [7]  from  the  calculus  of  residues  that 

.2 

Gl  -  £■[}£  [4cu,(())  cos  y  -  u"(0)]  (21) 


and 


TT  b  r  *  /rkv  1  2 

T  Uc  1“  (0)J 


(22) 


where  the  quantity  c  is  the  profile  chord  which  is  taken  as  unity  in  this 
work.  The  moment  can  be  calculated  after  the  complete  solution  has  been 
found . 


The  Form  of  m(£)  Near  the  Stagnation  Point 

This  is  also  a  well  known  result  which  we  shall  review  briefly.  The 
form  of  to  near  0  is  dominated  by  the  fact  that  on  a  smooth  contour  0  jumps 
by  the  amount  it.  In  particular,  as  one  traces  the  profile  surface, 
starting  at  Aj  in  Fig.  1  and  then  passes  through  0  while  proceeding  to  A£, 
the  jump  is  a  decrease  in  0.  This  is  precisely  the  behavior  exhibited  by 
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where  we  have  also  made  use  of  Eq.  (18)  because  it  applies  to  the  wetted 
surface.  But  from  Eq.  (15)  we  have 

T-in^-in/l-C  , 

U  p 

so  that 

.2  _ 

ds  -  2  -g—  [(cos  y  -  cos  B)  sin  8//1  -  Cp(B)]dB  (20) 

on  the  wetted  surface  of  the  profile.  As  we  have  noted  previously,  the  flow 
directions  differ  by  it  on  either  side  of  0.  Therefore,  if  the  sign  of  dz  is 
positive  on  the  arc  OA2,  it  will  be  negative  along  arc  0A^ .  As  a  result  of 
this  difference  ds  might  have  a  like  sign  change  in  these  two  regions.  Just 
how  this  might  occur  depends  on  the  form  of  /l  -  Cp(B)  in  any  particular 
case.  Therefore,  we  will  defer  consideration  of  this  question  to  a  later 
place  in  the  development. 

Hydrodynamic  Forces 

The  development  of  general  formulae  for  the  hydrodynamic  forces  on  the 
profile  depend  upon  certain  properties  of  the  function  co(c)  which  result 
from  the  previously  noted  fact  that  u>(c)  is  real  when  ?  is  real.  For 
then  one  can  apply  Schwartz's  principle  of  symmetry  in  order  to  write 
w(c)  =  <*1(5)  and  thereby  obtain  the  analytic  continuation  of  tu(c)  into  the 
lower  half  of  the  unit  circle  (6] .  Thus  we  can  write  for  a  prescribed 
modulus,  rQ  -  |?|  <  1, 

6(8)  -  ii(B)  -  fl(-B)  +  ix(-B)  , 
or 
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We  can  now  use  Eqs.  (13),  (14)  and  (15)  in  order  to  write 

.2  -r+ie  ,  . 

dz  -  dx  +  idy  =  — ^ -  [y  (?+“)“  cos  y](c  “  •  (17) 

On  the  wetted  surface  5  »  e1-^  and  Eq.  (17)  leads  to 

2b^  -t 

dx  =  e  [cos  y  -  cos  g]  sin  0  cos  8  dg 

and  ■  .  (18) 

2b2  -t 

dy  *  — —  e  [cos  y  -  cos  g]  sin  g  sin  0  dg 

Note  that  dy/dx  ■  tan  0  as  it  should  if  the  wetted  surface  is  to  be  a  stream¬ 
line.  On  the  upper  surface  of  the  cavity  t  *  0  and  arg  £  »  it  so  that  Eq.  (17) 
leads  to 

dxc  »  f-  cos  0  [\  (C  +  \) 

dyc  “  sin  9  [j  (?  +  f) 

provided  that  -  1  <  C  <  0. 

Returning  to  Eq.  (17),  we  can  write  the  square  of  the  arc  length  along 
the  wetted  surface  as 

2  —  0 
(ds)  ■  dz  dz  =*  {2  e  [cos  y  -  cos  g]  sin  g  dgj  , 


-  cos  Y](C  -  “•)  I2- 

(19) 

-  cos  y](c  -  £)  I2- 
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F  =■  b2[cos  Y  -  j  (? 


(in) 


Then  we  introduce  the  complex  velocity. 


,  dF 

w  =  u  -  iv  ■  —  , 

dz  ’ 


(ID 


in  order  to  write 


—  10  .2r  If  1  \  i  f  1  'i  1  d  £ 

w  -  qe  -  b  [cos  Y  -  y  (c  +  “ J  ](“  C  +  ^  ^7  » 


(12) 


or 


wdz  =  dF  =  b2  [j  (c  +  t*)  -  cos  y](?  “  ~  • 


(13) 


These  quantities  are  now  used  to  define  the  logarithmic  hodograph  or  (o-plane 

e_1“(?)  “  i  H  •  5  "  u  e"19  "  exp[£n  u  ~  i0J  *  (14) 


Therefore  we  have 


rn(^)  =  9  +  i£n  ^-  =  9  +  ix  ,  (15) 

where  t  =  In  On  the  free  streamlines  q  *  U  so  that  x  =  0  there.  In  the 
S-plane,  these  free  streamlines  are  on  the  real  axis  and  at  0*  we  know  that 
9=0  also.  Therefore, 

u(0)  =  0  (16) 

and  u)(s)  is  real  when  ?  is  real.  At  the  stagnation  point  q  ♦  0  so  that 
x  +  -  “  there.  The  flow  directions  differ  by  ir  on  either  side  of  0  and 
so  the  m-plane  with  the  various  corresponding  boundaries  can  be 
represented  as  illustrated  in  Fig.  6. 
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One  can  now  use  the  mapping 

/£7  -  /$, 

W  =  — -  (Z  +  cos  Y)  ,  (6) 

where 


cos  Y 


a  =  ^1  ~  /*2 
b  ^  $2  +  ^ 


(7) 


in  order  to  map  Aj  into  Z  =  +  1  and  A2  into  Z  -  -  1 .  The  point  0  maps  into 
Z  =  cos  Y  as  shown  in  Fig.  4.  The  significance  of  the  point  C,  also  shown 
in  this  figure,  is  discussed  later.  Then  the  Joukowski  transformation. 


z  -  -  J  U  +  ,  (8) 

maps  the  upper  half  of  the  Z-plane  onto  the  interior  of  the  unit  circle  in 
the  c-plane.  The  inverse  of  this  mapping  must  be 


c  -  - 


(9) 


in  order  to  make  the  point  0*  map  into  the  origin  of  the  C-plane.  Various 
corresponding  points  are  denoted  by  Aj ,  A2  and  0  in  Fig.  5.  The  arc  of 
the  unit  semi-circle  corresponds  to  the  wetted  surface.  Coordinates  on 
the  wetted  surface  are  given  by  ?  *  e*®  as  shown.  The  stagnation  point  is 
at  8  =  Y»  The  upper  surface  of  the  cavity  is  on  the  real  axis  between 
Aj  and  0*  and  the  lower  surface  of  the  cavity  is  on  the  real  axis  between 
A2  and  O'  as  marked  in  Fig.  5. 

One  can  now  use  Eqs.  (3)  through  (8)  in  order  to  write  the  composite 
of  the  preceding  mappings  as 
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sin  t, 


,2  sin  t,  cos  t 


yl  =  “ 


1  +  _T-L+  S1<ci) 


On  the  rest  of  the  wetted  surface,  that  is,  on  arc  OA2  between  the 
stagnation  point  and  the  tail,  we  know  that  8  <  y  and  we  put  t  **  -  t2» 
y  =  if  -  6  and  82  =  *  ~  ^2*  Then 


C2  b[cos  6  -  cos  ?2] 


In  this  case,  6  <  ^2  <  and  it  follows  that 


Then  we  have 


00  ^  f  <  - 

2  b( 1  +  cos  5) 


cos  t^  sin  t 


«(t_) 


sin  t„ 


y2  "  •“ 


+  -TT2t  si(t2)  -f 
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On  the  cavity  surface  AjO*  the  integration  starts  at  ;  =  -  1  and  ends 
somewhere  between  5  ■  -  1  and  C  »  0.  Note  that  the  value  of  t  at  5  *  -  1 
equals  that  of  t^  when  8  =  0.  So  that  if  we  put  t  =  t3,y  =  it-6  and 
C  =  -  C3,  we  have 


t 


3 


> 


then  the  upper  cavity  joins  the  nose  arc  of  the  wetted  surface  and  its 
coordinates  are 


£ 

U 


cos  t. 


sin  t„ 


+  Ci(t_) 


and 


(39) 


y3  =  — 


sin  t_  cos  t_ 

— T^  +  -T-1+  Si(t.) 


n 

2 


In  this  case 


0  <  t  <  rr-i - = - TT  • 

c  b(  1  -  cosT) 


Finally,  the  cavity  surface  from  the  trailing  edge  is  obtained  by 


putting  t*-t4,y  =  Tt-S  and  C  =  C4,  where  0  <  C4  <  1 .  Therefore 
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and 

0  <  t.  <  ,  , E - j-r  . 

4  b( 1  +  cos  6) 

Then  we  have 


and 


X4  “  U 


cos  t, 


y*  “  "  77- 


sin  t . 


+  d(tj 


sin  t. 


cos  t. 


+  Si(tJ  -  } 


(40) 


As  we  found  for  the  upper  cavity,  the  lower  cavity  surface  starts  smoothly 
from  the  wetted  surface. 

Equations  (37),  (38),  (39)  and  (40)  will  provide  the  shape  of  the 
wetted  surface  and  the  cavity  surfaces  for  the  simple  complementary  function 
of  Eq.  (33).  However,  these  equations  contain  the  undetermined  ratio  E^/U. 
We  shall  consider  the  ratio  E/b,  which  determines  the  strength  of  the 
complementary  solution,  as  a  parameter  which  we  can  prescribe  —  at  least 
for  the  time  being.  We  also  consider  the  value  of  y  (or  5)  to  be  known. 
Therefore,  we  need  to  "scale"  our  results  in  order  to  obtain  a  profile  of 
unit  chord.  Since  we  anticipate  that  the  complementary  function  can 
produce  a  rounding  of  the  wetted  surface  nose,  the  scaling  procedure  must 
account  for  this  possibility.  Accordingly,  we  shall  need  to  determine 
explicitly  the  location  of  the  apex  of  the  wetted  surface  nose  with  respect 


to  the  profile  chordline 
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Figure  1  shows  the  chordline  for  the  sharp-nose  or  for  the  round-nose 
case  when  the  upper  cavity  separates  at  the  leading  edge  with  respect  to 
the  hydrofoil  chord.  The  geometry  for  the  rounded  nose  with  the  separation 
point  on  the  upper  wetted  surface  behind  the  apex  of  the  wetted  surface 
contour  is  illustrated  in  Fig.  7.  The  apex  is  located  at  the  origin  of 
a,  n  coordinates  in  this  illustration  with  the  n  axis  being  tangent  to 
the  contour  at  this  point.  Denote  the  x,  y  coordinates  of  this  point  by 
xa»  Va*  Then  since  the  o-axis  is  normal  to  the  n  axis,  we  see  that  at  the 
apex  the  slope  of  the  contour  is 


y  a  rir  n 

■t—  =  tan  -  ctj  =  cot  a  . 


r\ 


We  will  restrict  our  attention  to  those  cases  in  which  the  apex  is  on  the 
arc  OAj .  Let  t  ■  tj  «  ta  at  the  apex.  Then  from  the  equations  preceding 
Eq.  (37),  which  define  tj ,  we  have 


'a  b(cos  £ -  cos 
a 


,  0  <  E  <  5  , 

u 


and  Eq.  (42)  can  then  be  used  in  Eqs.  (37)  to  define  xa  and  ya  once  the  value 
of  Ca  (°r  has  been  found.  Thus,  we  must  determine  the  unknowns  E^/U  and 
?a  in  terms  of  the  prescribed  quantities  E/b  and  <5(or  y) .  Two  conditions 
are  available  for  this  purpose.  The  first  is  given  by  Eq.  (41).  The  second 
will  be  that  the  profile  has  a  unit  chord. 

An  alternate  form  of  Eq.  (41)  is 

cos  a  =  sin  t  1 


sin  a 
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which  follows  from  the  complex  equation  just  above  Eq.  (36)  when  we  put 
t  »  ta.  We  will  use  Eq.  (43)  as  the  appropriate  expression  for  the  slope 
of  the  foil  contour  at  the  apex.  Now  let  us  differentiate  Eq.  (1)  g0 
that 


d(a  +  in)  =■  e^adz  . 

Then  from  the  complex  equation  just  preceding  Eq.  (36)  we  have 


e  iad(cr  +  in)  =  -  eiCdt  . 


(44) 


Starting  from  0  where  (a,n)  =  (a0,no),  we  integrate  Eq.  (44)  to  A2,  where 
(cr,n)  =  (1,0).  This  step  gives 


-ia 


0  -  «0  -  ‘"o)  -  r ' 


cos  t 


sin  t 

—  +  Ci('t) 


-  i 


sin  t  cos  t 


♦«(*.)  -t 


where 


Ct  b( 1  +  cos  6)  * 

corresponding  to  t2  with  ?2  "  11  in  Eqs .  (38).  Next  we  can  use  Eqs.  (37)  with 
?1  *  ?a  and  tj  -  ta  or 


JCuJt.  -  ■»  ’-*  *-*  *-•  •  ‘ 
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£  s  .  ,  ■  .  ,  ■  .  ■  ■ 

a  b(cos  £  -  cos  6) 

a 


at  the  apex  in  order  to  integrate  Eq.  (44)  from  0  to  the  apex  along  the  arc 
OAj .  This  step  results  in 


-  e 


-ia 


K  +  iT1o^  =  U“ 


cos  t  sin  t 
a  a 


„  +  Ci(t  ) 

t  '  a' 


+  i 


sin  t  cos  t 

~~r£’  * 

t  a 


TT 

2 


Eliminating  the  sum  Oq  +  iriQ  from  these  two  equations,  we  get 


and 


cos  a  -  F  -  f(t  ) 


—■  sin  a  =  G  -  g(t  ) 
El  3 


(45)’ 


In  Eqs.  (45)  we  have 


*The  F  introduced  here  is  not  to  be  confused  with  F  =  <j>  +  ii|>  from  Eq.  (2). 
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and 


F 


xtU  cos  Ct  sin 


sin  tt 


2 

t 

t 


COS  tt 


t 

t 


(46) 


which  contain  known  quantities  because  tt  is  known.  The  remaining  pair  of 
functions. 


cos  t  sin  t 


a 


'i 


and 


(47) 


sin  t 

a 


+ 


cos  t 

a 


t 

a 


contains  ta  which  depends  upon  the  unknown,  Ba  *  tt  -  Ca.  Thus,  Eqs.  (45)  are 
two  simultaneous  equations  containing  the  quantities  U/E^,  a  and  8a  which 
must  be  determined.  Therefore,  Eqs.  (43)  and  (45)  form  a  determinate  system 
which  can  be  solved  by  iteration.  In  order  to  do  this,  we  can  write  the 
complete  system  as 

G  "  g(ta) 

F  -  f[t  J  ’ 

^  a' 


tan  a 


(48) 
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cot  t  =  tan  a 
a 


(49) 


and 


F  -  f([,) 

cos  a 


G  -  8^a) 

sin  a 


(50) 


In  the  derivation  of  this  system  we  have  assumed  that  the  apex,  za,  is^ 
on  the  arc  OAj .  On  the  other  hand,  we  specify  the  quantities  E/b  and 
y  *  tt  -  S .  We  must  now  determine  whether  or  not  our  assumption  regarding 
the  location  of  za  can  restrict  possible  choices  for  the  parameters  E/b  and 
6.  In  particular,  we  recall  that,  as  is  true  for  the  quantity  Cj,  we  must 
also  require  that  0  <  £a  <  5,  as  noted  in  Eq.  (44).  The  limiting  condition, 
corresponding  to  the  coincidence  of  the  apex  and  cavity  separation  point  at 
the  nose  of  the  profile  corresponds  to  Ca  =»  0.  In  this  case,  the  smallest 
value  of  tj  for  any  choice  of  E/b  will  be  found  when 


t 

a 


min 


E _ 

b(l  -  cos  6) 


On  the  other  hand,  by  inspection  Eqs.  (46)  and  (47),  we  see  that  the  largest 
values  of  f  and  g  are  found  for  ta  =  tm£n.  The  values  of  F  and  G  are  also 
obtained  from  the  smallest  value  of  t2  because  tt  is  calculated  from  t£ 
with  ^2  =  >  namely: 

t  _  E 

t  b( 1  +  cos  6 ) 

Let  us  compare  the  values  of  F  with  f  and  G  with  g.  Suppose  that  E/b  is 
selected  so  that  cos  ta  "  cos  tt  *  1 ,  sin  ta  *  ta  and  sin  tt  B  tt.  Let 
5  <<  1 .  In  this  case  we  can  see  that 
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* 

t 

a 

•* 

min 

II 

C"4 

•  \ 

1  - 

and  —  - 

■ 

t 

a 

min 

-  f-)2 

,-ii 

Ct 

1  3 

ana  — 

g 

Ct 

( 

k 

l 

. 

For  example,  if  6  =  .1,  we  would  estimate  F/f  ~  (20)^  and  G/g  ~  (20)2.  These 
estimates  imply  that  both  E/b  and  6  are  significantly  smaller  than  unity. 

In  the  applications  contemplated,  E/b  will  probably  be  less  than  unity 
although  6  might  conceivably  approach  or  exceed  unity.  Therefore,  we  shall 
consider  the  ratio, 


t 

a 


min 


t 

t 


n 


1  +  cos  6 
1  -  cos  6 


n  >  1  , 


which  permits  us  to  consider  roughly  the  ratios  of  F  to  f  and  of  G  to  g 
for  various  values  of  6.  In  particular,  we  can  solve  for  cos  6  and  obtain 


cos 


6 


n  -  1 
n  +  1 


> 


which  permits  us  to  plot  a  curve  of  n  =  (t  /t  )  versus  6  as  shown  in 

3 'min  C 

Fig.  8.  This  curve  illustrates  the  effect  that  the  choice  of  stagnation 
point  location  has  on  the  ratio,  n.  The  value  of  n  in  turn  gives  a  rough 
indication  of  how  large  the  ratios  F/f  and  G/g  will  be. 

It  appears  for  most  cases  that  these  ratios  will  be  very  large  and 
one  need  not  solve  Eqs.  (48)  and  (49)  by  iteration.  Instead  one  can 
obtain  an  accurate  value  of  a  from 


►  t.  G 

tan  a  *■  — 


(48a) 
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and  he  can  then  determine  U/E^  from 


U_  _  F 
2  cos  o 

£ 


(50a) 


Should  cases  arise  in  which  Eqs.  (48a)  and  (50a)  are  not  accurate,  they  can 

be  used  advantageously  to  start  the  iteration.  In  order  to  illustrate  these 

points  and  in  order  to  show  a  profile  shape  derived  from  the  complementary 

solution,  we  have  prepared  the  following  numerical  example. 

We  started  the  calculation  by  selecting  6  =  70°  and  E/b  =  .01. 

Figure  8  shows  that  n  *  2.  The  values  t  =  .01520  and  t  =  .00745 

a  ,  t 

min 

follow  from  the  formulae  for  these  quantities.  From  Eqs.  (46)  we  find  that 
F  =  18,004  and  G  ■  266.83.  Equations  (48a)  and  (50a)  lead  to  tan  a  =  .01482 
and  U/E^  =  18,006.  From  the  formulae  just  after  Eq.  (35)  we  have  for  the 
cavity  drag  due  to  a  profile  of  unit  chord 


CD  =  2ir  =  .00035 
c 


as  the  contribution  for  this  point-drag  profile.  The  value  of  ta  can  now 
be  found  from  Eq.  (49).  It  is  ta  =*  1.556.  Equation  (42)  can  now  be  used 
to  find  that  5a  =  69.61°.  Note  that  for  this  case  the  apex  is  almost 
coincident  with  the  stagnation  point.  Because  Eq .  (34a)  shows  the 
complementary  function  to  be  at  the  stagnation  point  this  result  is 
expected.  The  fact  that  the  apex  is  not  exactly  at  5  =  70°  is  due  only 
to  the  inclination,  a  =  .85°,  between  the  chord  line  and  the  x  axis. 

Once  U/E^  has  been  found,  the  values  of  xt  and  yt  can  be  found  from 
Eqs.  (46).  The  values  are  xt  =*  .99989  and  yt  =*  -  .01482.  Now  one  can 
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use  the  conditions  that  in  Eqs.  (1),  (xt,yt)  <-►  (a=l,n=0)  and 
(xa,ya)  «-*■  (o=0,n=0) .  These  lead  to  a  system  of  equations  from  which 
0g  and  ng  can  be  eliminated  and  one  finds  that 

x  =  x  -  cos  o  =  -  sin  t  , 

at  t  a  * 

y_  =  y.  +  sin  a  =  y_  +  cos  t 

J  a.  J  t  J  t  a 

When  the  above  values  of  xt  and  a  were  used  in  these  equations,  the  values  of 
xa  and  ya  were  found  to  be  zero  to  within  five  decimal  places.  This  result 
is  consistent  with  the  location  of  Ca  noted  previously.  Continuing  with 
Eqs.  (1),  we  can  use  the  fact  that  xa  =  ya  =  0  to  see  that  it  must  also 
follow  that  <*0  33  ng  33  0.  Accordingly,  the  form  of  Eqs.  (1)  for  the  present 
calculations  is 

o  =»  .99989x  -  .014796y  , 
n  33  .014796x  +  .999989y  . 

The  next  phase  of  the  calculations  is  the  evaluation  of  the  equations 
for  the  wetted  surface  and  cavity  contours  in  accordance  with  Eqs.  (37), 

(38)  and  (39).  The  result  of  these  calculations  is  shown  in  Fig.  9.  In 
this  figure,  the  chordline  distance,  a,  has  been  labeled  as  X  and  the 
ordinate,  n,  has  been  labeled  as  Y.  Note  that  the  Y-scale  is  magnified 
five  times  compared  to  that  of  the  X-scale.  The  trailing  edge  of  the 
wetted  surface  is  at  X  3  1.  The  upper  surface  separation  point  is  at 
a  »  X  =  .240.  The  cavity  thickness  at  X  3  1  is  Y  =  T  3  .02980.  This 

point  is  marked  to  the  same  scale  as  the  X-scale  by  the  dot  and  the  line 

at  X  3  1  in  order  to  give  an  idea  of  the  actual  thickness  of  this  example 

of  a  point-drag  profile.  Finally  we  can  calculate  the  value  of  ?c(l) 
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in  this  example  in  order  to  compare  it  with  the  values  found  previously  for 
the  flat-plate.  It  is  found  from  the  value  of  tc  =  (E/b)/[(£c  +  l/^c)/2  cos  y 
at  cr  =  1 .  The  calculations  indicated  give  Cc(l)  =  .3290  which  is  roughly  two 
times  the  values  of  Table  1. 

Our  desire  to  retain  as  much  simplicity  as  possible  in  the  preceding 
analysis  of  the  complementary  function  has  caused  us  to  place  the  point-drag 
singularity  at  the  stagnation  point.  This  restriction  on  the  location  of  the 
point  of  application  has  allowed  us  to  show  that  such  a  solution  exists, 
that  it  definitely  leads  to  a  smoothly  rounded  profile  nose  and  that  it  will 
cause  an  incremental  thickening  of  the  cavity  depending  on  its  strength,  E. 

Of  course,  we  need  not  restrict  ourselves  to  the  stagnation  point  as  being 
the  location  of  the  point-drag  singularity. 

For  example,  suppose  we  choose  some  other  point  C  on  the  wetted  surface. 
Such  a  point  is  illustrated  in  Fig.  5  and  it  happens  to  be  located  between 
the  upper  cavity  separation  point  and  the  stagntion  point,  although  C  could 
just  as  well  be  at  some  other  wetted-surf ace  location.  The  main  idea  is 
that  now  6  =  yc  at  the  location  of  the  point-drag  signularity  and  if  we 
simply  replace  Eq.  (34a)  by  the  modified  expression, 


“c<C> 


-  2E t. 


b[cos  Yc  -y(c  +  £)]  b(c  .  eiYC)(i;  _ 


-iy. 


(34b) 


we  still  have  a  function  which  satisfies  those  conditions  needed  for  a 

complementary  solution.  It  is  clear  that  in  the  £-plane,  a>c( 1 )  =  0. 

Moreover,  when  5  is  real  ojc  is  real  and  on  the  unit  circle  rc  =  0  every- 

+  iy 


where  except  possibly  at  the  simple  poles,  5  =  e 


From  Eqs .  (21)  and 


(22)  it  follows  that 
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(51) 


as  before.  On  the  other  hand,  because  of  the  displacement  of  the  point  C 
away  from  0,  a  lift  force  is  produced  and  we  find  that  now 


8nbE 

Uc 


Y  +  Y, 


sin 


sm 


8TTbE 

Uc 


sin 


6+6 

c 

2 


sin 


,  (52) 


where  Y  ■  it  -  6  and  Yc  =  accordance  with  previous  convention.  The 

profile  chord,  c,  should  be  set  at  unity  in  Eqs.  (51)  and  (52).  Equation  (52) 
shows  that  CL  =  0  when  6C  =  6.  However,  6C  <  6  when  the  point  C  moves  toward 
the  point  Aj  and  a  negative  lift  results.  In  the  limit  as  6C  +  0,  we  have 


8TfbE 

Uc 


sin 


2  6 
2 


If  C  is  between  0  and  A2  a  positive  lift  is  produced  and  in  the  limit 
when  C  is  coincident  with  A2  we  have 


8irbE 

Uc 


2  6 

cos  j  . 


If  one  were  to  let  E  be  negative  the  sign  of  the  foregoing  trends  with 

respect  to  CL  would  be  reversed.  However,  we  must  insist  that  E  >  0 

because  this  function  produces  a  thickening  of  the  cavity  and  because  then 

a)  (0)  =  -  2E/b.  Similarly,  we  have  also  found  that  wn(0)  =02  sin  y. 
c  u 

Thus,  the  effect  of  adding  u)q  and  u>c  increases  the  net  drag.  Neither  of 
these  functions  can  act  to  reduce  it.  Accordingly,  we  shall  take  Eq.  (34b) 
as  the  appropriate  form  of  the  eigensolution  which  has  been  sought . 


-36- 


25  January  1985 
BRPrlhz 


Conclusions 

The  chief  finding  of  this  paper  is  that  one  can  construct  many  eigen- 
solutions  for  the  exact  inverse  problem  of  two-dimensional  cavity  flow  at 
zero  cavitation  number.  From  among  these,  we  have  chosen  that  eigensolution 
which  appears  to  offer  the  least  disturbance  to  the  upstream  flow  field 
inclination  of  any  cavity  flow  which  does  not  already  include  a  point-drag 
solution  as  one  of  its  elements.  This  particular  choice  also  seems  to  offer 
the  greatest  analytical  convenience.  The  physical  conditions  satisfied  by 
this  eigensolution  are: 

(1)  At  points  on  the  cavity  and  on  the  wetted  surface  of  the  profile, 
the  flow  velocity  is  equal  in  magnitude  to  the  free-stream 
velocity . 

(2)  The  point-drag  solution  vanishes  at  infinity,  but  it  does  have 
an  integrable  singularity  on  the  profile  surface. 

(3)  This  function  produces  no  additional  flow  inclination  on  the 
entire  upstream  stagnation  streamline. 

A  specific  example  of  the  flow  geometry  represented  by  an  isolated 
eigensolution  has  been  given  above  to  show  how  this  function  can  produce 
round-nosed  profiles.  In  general,  it  is  found  that  the  point-drag  solution 
produces  a  widening  of  the  cavity  which  is  directly  proportional  to  its 
strength.  An  incremental  cavity  drag  accompanies  this  widening  and  this 
drag  is  proportional  to  the  square  of  the  eigensolution  strength.  No  lift 
is  produced  by  the  point-drag  function  when  its  location  coincides  with 
that  of  the  stagnation  point  on  the  profile  surface.  In  contrast  to  the 
linearized  theory,  the  complementary  function  singularity  need  not  be  at 
the  stagnation  point.  In  these  cases,  the  incremental  cavity  drag  is 
not  changed  from  its  value  when  the  singularity  is  at  the  stagnation  point. 
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But  when  the  singularity  is  located  between  the  stagnation  point  and  the 
upper  separation  point  a  negative  incremental  lift  is  produced.  If  the 
singularity  is  on  the  lower  surface,  downstream  of  the  stagnation  point, 
a  positive  lift  increment  is  found. 

As  a  result  of  these  findings,  it  appears  that  an  eigensolution  exists 
for  the  nonlinear  theory  of  cavity  flow  at  zero  cavitation  number  and  that 
it  is  now  most  liekly  that  a  similar  eigensolution  can  be  found  for  such 
fully  cavitating  flows  at  cavitation  numbers  which  are  greater  than  zero. 
The  results  found  so  far  suggest  that  the  nonlinearized  theory  and  the 
linearized  theories  parallel  one  another  very  closely  as  far  as  the  nature 
of  the  point  drag  solution  are  concerned.  But  the  present  results  exhibit 
some  features  which  are  lost  whe.i  the  process  of  linearization  is  applied 
in  the  linearized  formulation.  Finally,  the  way  is  now  clear  to  augment 
the  method  of  Levi  Civita  by  adding  in  the  eigensolution,  wc,  in  order  to 
formulate  an  exact  inverse  theory  which  permits  one  to  control  the  cavity 
thickness  and  to  prescribe  the  wetted-surface  pressure  distribution 
simultaneously. 
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Cameron  Station 
Alexandria,  VA  22314 
(Copies  1  through  6) 

Commanding  Officer 
David  W.  Taylor  Naval  Ship 
Research  &  Development  Ctr. 
Department  of  the  Navy 
Bethesda,  MD  20084 
Attention:  M.  Tod  Hinkel 
Code  1504 
(Copies  7  through  12) 

Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 
Attention:  T.  E.  Peirce 

Code  NSEA-63R3 1 
(Copy  No.  13) 

Director 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16804 
Attention:  M.  L.  Billet 

(Copy  No.  14) 

Director 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16804 
Attention:  L.  R.  Hettche 
(Copy  No.  15) 

Director 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16804 
Attention:  J.  W.  Holl 
(Copy  No.  16) 


Director 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16804 
Attention:  B.  R.  Parkin 

(Copy  No.  17) 

Director 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16804 
Attention:  GTWT  Files 
(Copy  No.  18) 

Director 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16804 
Attention:  ARL/PSU  Library 
(Copy  No.  19) 
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